Vector Nonlinear Klein-Gordon Lattices: General Derivation of Small Amplitude 
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Group velocity and group velocity dispersion for a wave packet in vectorial discrete Klein-Gordon 
models are obtained by an expansion, based on perturbation theory, of the linear system giving the 
dispersion relation and the normal modes. We show how to map this expansion on the Multiple 
Scale Expansion in the real space and how to find Non Linear Schrodinger small amplitude solutions 
when a nonlinear one site potential balances the group velocity dispersion effect. 

I. INTRODUCTION. 



One of the most popular approaches used to determine the small amplitude envelope soliton solutions in non linear 
models is the well-known Multiple Scale Expansion (MSE) technique This technique amounts to expanding the 
equations of motion on different time and space scales looking for wave-packet-like solutions; a wave packet is a 
superposition of plane waves whose frequencies and wave vectors lie in a narrow band, and it can be conveniently 
CO , described by a plane wave with an amplitude that varies slowly in space and time. Increasing progressively the 
■ time and space scales one determines in a first step the carrier wave as a phonon mode of the linearized system, 
' then deduces the partial differential equation that identifies the envelope phase velocity with the wave packet group 
velocity, and finally derives the NLS equation for the envelope whose diffusion coefficient is in fact the wave packet 
I group velocity dispersion. An alternative method, commonly used in optics, consists in expanding the dispersion 
relations with respect to the carrier frequency and then in building at each order of this expansion an operator that 
acts on the envelope function ||^. MSE has been successfully applied to various non linear systems with scalar fields 
and the corresponding method adopted in optics has permitted the study of optical solitons in fiber. In this latter case 



I one deals with eletric field components which are not coupled at the linear order of the Maxwell equations. However, 
O many non linear models of interest involve vectorial fields with coupled components at the linear order that give rise 
^ to dispersion relations with more than one branch: classical examples are given by multi-atomic lattices, or by lattices 
'■^ in which the mass at each site can move in a multidimensional space. 
• • In this work, we show how to find small amplitude envelope soliton solutions in such vectorial lattices problems. 

The main difficulty with respect to the scalar case is to determine the relative amplitudes of the different components 
of the field. We will perform a perturbative expansion, around one wavenumber, of the linear system that gives the 
^ ' dispersion relations and the linear eigenmodes; then we will introduce an operator formalism analogous to that used 
5^ in optics to obtain, from this expansion and from the nonlinear terms, the MSE equations, up to the NLS one. 

An application of the method presented in this work can be found in a forthcoming paper where the envelope 
soliton solutions of an hclicoidal DNA model described by a radial and an angular degree of freedom for each site are 
derived. 

II. WAVE-PACKET IN LINEAR VECTORIAL LATTICES 

The NLS equation is obtained when a weak dispersion is balanced by a weak nonlinearity. In order to characterize 
the dispersion, let us first restrict our attention to the linear part of the system of interest. We start with a one 
dimensional vectorial linear lattice model given by the equations of motion: 

= - J{n-n ,a,a)E{n ,a ,t) (1) 

n' ,a' 

where n,n' are the site indices, a, a' are the indices that label the components of the vectorial field E{n,a,t), and 
J{n — n' , a, a') are the force constants depending on n — n' for translationally invariant systems. 
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Looking for plane wave solutions of the form 

Al^(g)e*(«"-'^'(«)*) +C.C. (2) 

where A is the wave amplitude, the equation of motion is mapped to the operator equation in the wave numbers 
space: 

{J{q)-u;f{q))Vi{q)=0 (3) 

where J{q) is the Fourier transform of the matrix J{n — n'). The index / runs from 1 to the number of components of 
the vectorial field E{n,a,t)] the eigenvalues functions Luf{q) give the branches of the dispersion relation; the normal 
modes Vi{q) are the orthonormalized eigenvectors of the matrix J{q) — ujfiji). 

In order to investigate the dispersion, we now consider a wave-packet- like solution, i.e. a superposition of plane 
waves with wave numbers in a small interval: 

El (n, t)^ A{q) Vi (q) e*(«"-'^' dq + c.c. (4) 

For each q contributing to the wave packet the system of equations (||) must be fulfilled. The weakly dispersive 
case is obtained by considering only small deviations of q with respect to the wavector qg corresponding to the center 
of the wavepacket. To measure this deviation, the wavevector q is written q — qo + £qi, where e 1. 

Eq. is solved, Vqi in the integration range, by a perturbative technique. 

The operator J{qo + eqi) can be expanded in Taylor series as J{qo) + eJ'((;o)9i + J" {<io)<l\/'^ + . . .. The quantities 
eJ'{qo)qi, e^J"{qo)ql/2 are small perturbations with respect to unperturbed operator J{qo), whose eigenvalues are 
ujf — iLifiqo) and whose eigenvectors Vi — VJ (go) constitute a complete basis. 

According to standard perturbation theory [Bj we write the expansions of the eigenvectors and eigenvalues, 

Viiqo + eqi) = Vi + eV.^'^qi + e%'^^'>ql/2 + ... (5) 
^liqo + egi) =L^i+ + e^ujPql/2 + ... (6) 

Equation (^) has to be solved at each expansion order: 

at order e° : JVi = ujfVi (7) 
at order : (J\^^^^ + J'Vi) qi = {2u;iuj'i'^^Vi + V^^^^) qi (8) 

at order : (J'V,^'^ + \j"Vi + \jV^^^) ql = (c^^^'y, + -I- 2c.,c.Wy/^) -t- ^V^^^^) qf . (9) 



At order e*^, one solves the unperturbed problem determining Vi and uJi. At order e, one determines Vi^\ uj[^^ 
(1) 



UJ 



imposing to to be orthogonal to Vi to guarantee the normalization of V;(go + e^i) the scalar product of (g|) with 
^ gives 

vl"-^ = o^^y^ (10) 
y^j'Vi 

am = ^ T (11) 



and that with V* gives 



At order e^, one determines cjp'' by multiplying (^) by VJ* 



UJ 
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We assume, for sake of simplicity, that the eigenmodes of J (see (||)) are not degenerate, but the generaUzation of the 
degenerate case is straightforward with the standard perturbation theory. 

The phase of each component of (||) can be expanded around the central wave number qq, up to second order in 

eqi = q — qo using the values of wp"* and ojp^ determined above 
^z(n,t)=e^(«°"-'^'(''°)*)e / Aiqo + eq,)Vi{qo + eq,) exp \ieq,{n - iul'\qo)t) - ie^^iol^\qo)t\ dq, + c.c. (14) 



Under this form, Ei(n, t) appears as a plane wave, henceforth called the carrier wave, with an amplitude that depends 
on space and time and which corresponds to the integral of equation (^^, Ei{n, t) — F{n, t) exp[i{qQn — uji{qQ)t)]+ c.c. 
The fact that ijj\^\qf)), ijj'f'\qQ) obey relations ([l^ ) and ( [l^ ensures that this wave packet is a solution of the original 
equation (^, up to the order of the various expansions. 

In order to extend the study to the nonlinear case it is useful to express these conditions under the form of an 
equation in the space time coordinates for the amplitude. Let us introduce the quantity 

/+A9/e r 2 > 

^ ^ A{qo + tqi) exp|iegi(n-w^')(go)t) - y^r(go)ij dq^ (15) 

Equation ( |l5| ) shows that A{n^ t) slowly varies in space and time. In the spirit of the multiple scale expansion, it is 
natural to introduce the slow variables Xi — en, ti = et and ^2 = so that A(n, t) can be written as 

A(go + egi) exp|igi(a;i - W;*^^(go)ti) - «yt^P' ('7o)i2 1 dqi (16) 



A{qQ + eqi) exp jigisi - iywp^ (90)^2 1 dqi (17) 



with the introduction of the variable si — xi — cj^'' ((70)^1 to switch to the frame moving at the group velocity of the 



carrier wave 



Using the relation {dA/dxi) — {dA/dsi) — i < qi >= J^^qj^ iQi A{qQ + eqi) exp |?(?iSi — i^i^P^ (90)^2 1 dqi that 

derives directly from Eqs. ([l6| ) and ([l7|), and the expansion (||) of V{qo + eqi), the amplitude F of the wave can be 
expressed as a function of A{si, 12) by the relation 



F{xiMM) = <Vi - ieV{^)A{xi,ti,h) . (18) 

0x1 



From (|16D and ([I3), we directly derive the equations of motions of ^ as a function of the slow space-time variables: 



and 



.OA d^A 



dt2 2 dsi 



= (20) 



where and wp-* are then the group velocity and the group velocity dispersion of the wave packet and determine 

the peak velocity and the spread out of the envelope function. Equation (|lj) shows that V'l determines the first order 
correction to the direction of the vectorial field solution. 



III. NON LINEAR VECTORIAL LATTICE 



We now consider the full equation of motion, including nonlinear on-site potential terms. Extra nonlinear terms 
depending on time derivatives may also appear in the case of non Cartesian coordinates systems: 
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0^ = — J{n - n ,a,a )h{n ,a ,t) + 

n' ,a' 

d—0 k—Q a" ,a' <.a" 

E E E E ^x^-.. ""')^'-''^ ^)^''"-''^ ('^' ^)^''"'^ ^) (21) 

/c— J — a'" ,a" <a"\a' <a" 

where E'^^\n,a^t) indicates the j—th time derivative of i?(n, a, i), and cj| j,(a', a"), i. j{a' ,a" ,a"') are the quadratic 
and cubic nonhnear terms numerical coefhcients. Index d is the total time derivative order of each term. The terms 
with d — are the nonlinear potential force terms; the others derive from the kinetic energy in the case of non 
Cartesian coordinates so that d < 2. 

The quadratic terms in ( pl| ) give rise to second harmonic and constant terms that, if we look for small amplitude 
solution, have to be included as additional smaller corrections: 

+ e2e2'('°"°-"'*")7/^'(a;i,ii,t2) + \A{x,,h,h)\^ + Oie"") . (22) 

We are interested in situations where dispersion can balance nonlinearity, and therefore have to be measured by the 
same scaling parameter e. While the overall e factor was not important in the linear case, it must be explicitly included 
in the nonlinear case. 

We solve the equation of motion ( pT| ) on the three characteristic magnitude scales of the wave packet. At order e 
we get ui ,Vi from equation (^. 

At order we get for the wave packet term expressed in the form the system of equations (|^) for each qi in the 
integration range. After integration on the envelope distribution this gives rise to the equation in the anti-transformed 
Fourier space 

(2..^ A + (i - .f)^/^) A + Ai.,M) . . (23) 

In fact from ([l^), the average wave numbers deviation is < qi >= —idA/dxi and the averaged frequency deviation 

is in the same way < Alui >=< ujI^\i >= idA/dti. 

By scalar product of (23) and Vm ^ I we obtain the components a„i ([ll] ) of V^^'^'' ( |l0| ) on the base {Vm} and by 

scalar product of ( p3| ) and V^* wc obtain the equation (|l9|) with W;"'^'' defined by (p^). 

At the same order of expansion one determines the vectors 7;, fli collecting the terms of corresponding order (e^) 
and phase in the equation of motion ( |2l| ) in which the solution form (j2^) has been inserted [^. One then obtains 7; 
by solving the algebraic system 

2 d 

(J(2go)-W(go))7i=EE E Cd,fc(a',a")(-*c.z)'^(«')Vi(a") (24) 

d=0 k=0 a" ,a'<a" 

where Cc;,fc(a', a") is the vector of components ^^(a', a"), each derivative with respect to to giving a factor (— iw;). 
And \ii is obtained from the system 

2 d 

'hmi-^Y^Y. E Qi,fc(a^a'')[(^c.^)'(-*^^)'"'^*("')^(«'') + (-*^^)'(*^^)'''^^*(«'')^^(«')]• (25) 

rf— fc— a." ^ol' <ol" 

At order e"^, from (|2^) and ( p2| ) for the terms in e*^'^"""~^^^'^"^*"\ one obtains the system of equations: 

'[(i J" - c.p)')l^ + (J' - 2c.,c.«)l^(i) + (^(-^ - )^^'')]^ + 2*^'^'^) ^(^1' ^2) + 

Q|A(si,t2)P^(si,t2) = (26) 

where 
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2 = EE E Q,fc(a^a")[(^c.0'(-2^^^)'"'^;(«')7^(«") + {-2iu;i)\tu;if-'^V: {a")^i{a')] + 

d—0 k—0 a" ,a' <.a" 

2 

2E E Cd,o{c'',a"){-iujiffii{a')Vi{a") + 

d—0 a." .a' <.a." 

2 d k 

EEE E C,,,,,(a',a",a"')[(-*^;y(-*^0'"-'(*^0'~'^K«')^(«")^*K") + 

d=0 fc— i— a'" ,a" <.a'" .a' <.a" 

{-iiUi)^iu;i)''~^{-zu;if-''Viia')Vria'')Viia''') + i^^^^^^ (27) 

The first line of the equation ( p^ ) corresponds to the third order expansion of the linear operator equation ^ 
applied to the wave packet (terms in e''i'iono-^i{qo)to) ([22[)), in the moving reference frame used in (pO|), < qi >= 

-i{d/dsi)A{si,t2), and < ioPqf/2 >-< A(A(a;()) >- i(a/9i2)A(si, ia). 

The first two lines of the non linear coefficients vector Q arise from the double product, in the nonlinear quadratic 
force terms in (pl| ), between 0(e) and O(e^) components of (p^); the last two lines come from the nonlinear cubic 
force terms in ( |2l| ) when considering just the 0{e) terms in Ei{n,t). 

Multiplying Eq. (|2|) by V^* we obtain the NLS equation: 

.d_ 

where 



^7J:2 +«— ) MsiM) + Q\A{si,t2)\^ A{siM) = (28) 




and 



V*Q 

The first part of equation (^) is equation (^^ for the wave packet with 2P = cjj^-* given by the (^3|) to which is 
now added the nonlinear part with coefficient Q. If PQ > then the effect of the amplitude dependent non linear 
potential well in ( |28| ) balances the wave packet group velocity dispersion giving rise to the stable envelope soliton 
solution 0]. 



IV. SUMMARY. 



The main outlines of the approach we introduce are the following. The envelope soliton like solutions arise in systems 
with weak dispersion and weak nonlinearity by two parallel series expansions driven by a common expansion parameter 
(e): on one hand the weakness of the diffusion, for a wave packet like solution, allows an expansion of the equations 
that regulate the space time behaviour of the solution on different scales; on the other hand the weak nonlinearity, 
for small amplitude solutions, allows to write the equations of motion at increasing orders of accuracy introducing 
the nonlinear terms in a progressive way. For scalar fields the Taylor series expansion of dispersion relations gives 
directly the diffusive part for the envelope equations of motion. Vectorial fields are instead characterized by a linear 
part which gives rise, in the q space, to an eigenvalue (dispersion relations) and eigenvector (relative amplitude of the 
different components) problem (|^) ; to obtain the correct expansion in multiple scales it is then necessary to apply the 
perturbation theory (^,(^). Finally, to combine this perturbative expansion with the nonlinear one, we antitransform 
the equations (||), (|9|) as done in (p^), (p6|). 

Following the approach introduced in this paper it is straightforward to derive the NLS equation for every non 
linear vectorial lattice with on-site nonlinearities and with an arbitrary number of components. For more complex 
systems this could even be programmed in symbolic languages to provide a fully automatic method. After having 
identified the nonlinear coefficients ,a") and C2 ^ ,a" ,a"') in the equation of motion (pl|), there are only 

algebraic systems to solve: one has to solve the eigenvectors Vi(qq) and the eigenvalues uji{qo) of the matrix J{qo), 
then the systems ( p^ ) and (^5|) for 7; and , and to derive P from (|2^) and Q from (^^ and ( |30| ) . From (|2|) one then 
obtains, if PQ > 0, the envelope function A{xi,ti,t2) that, inserted into ( p^ ) together with the eigenvector correction 
(p^, (pi]), gives the complete O(e^) solution we are looking for. 
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